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Abstract: The double cone, a cone over a product of a pair of spheres, is known to play 
a role in the black-hole black-string phase diagram, and like all cones it is continuously self 
similar (CSS). Its zero modes spectrum (in a certain sector) is determined in detail, and it 
implies that the double cone is a co-dimension 1 attractor in the space of those perturbations 
which are smooth at the tip. This is interpreted as strong evidence for the double cone being 
the critical merger solution. For the non-symmetry-breaking perturbations we proceed to 
perform a fully non-linear analysis of the dynamical system. The scaling symmetry is used 
to reduce the dynamical system from a 3d phase space to 2d, and obtain the qualitative 
form of the phase space, including a non-perturbative confirmation of the existence of the 
"smoothed cone". 
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1. Introduction and Summary 

The phase diagram of the black- hole black-string transition (see the reviews Q) was 
conjectured in Q to include a "merger" point - a static vacuum metric (see figure ^ which 
lies on the boundary between the black-string and black-hole branches. It can be thought 
as either a black string whose waist has become so thin that it has marginally pinched, or 
as a black-hole which has become so large that its poles had marginally intersected and 
merged (and hence the name "merger"). The metric cannot be completely smooth as it 
interpolates between two different space-time topologies, but it may have only one naked 
singularity. 

The black-hole black-string system has been the subject of intensive numerical research 
Hi ^! Hj 05 H • Naturally, the merger space-time itself is unattainable numerically since it 
includes a singularity, but it may be approached by following either of the two branches far 
enough. Indeed, all the available data indicates that the black-string and the black-hole 
branches approach each other, in accord with the merger prediction. 

At merger the curvature is unbounded around the pinch point. One defines the "critical 
merger solution" to be the local metric around the pinch point (at merger), namely the 
one achieved through a zooming limit around the point. It is natural to predict 1^, ^ that 
the critical merger solution will lose all memory of the macroscopical scales of the problem 
(the size of the extra dimension and the size of the black hole) and moreover be self-similar, 
namely invariant under a scaling transformation. The central motivation of this paper is 
to determine the critical merger solution. 

Self-similar metrics belong to one of two classes: Continuous Self-Similarity (CSS) or 
Discrete Self-Similarity (DSS): while a CSS metric is invariant under any scale transforma- 
tion and can be pictured as a cone, a DSS metric is invariant only under a specific scaling 
transformation (and its powers) and can be pictured as a wiggly cone with its wiggles being 
log-periodic (see figure §). A key question is: Is the critical merger solution CSS or DSS? 
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Figure 1: The merger metric, r is the radial coordinate in the extended directions, z is periodically 
compactified, while time and angular coordinates are suppressed. The heavy lines denote the horizon 
of a static black object which is at threshold between being a black-hole and being a black-string. 
The naked singularity is at the x -shaped pinching (horizon crossing) point. Upon zooming onto the 
encircled singularity it is convenient to replace (r, z) by radial coordinates (p, x) radial coordinates. 
We shall be mostly interested in the "critical merger solution" - the local metric near the singularity, 
namely the encircled portion of the metric (in the limit that the circle's size is infinitesimal). 



The significance of the critical solution 
is that critical exponents of the system near 
the merger point are determined by properties 
of this solution (this is known to be the case 
in the closely related system P] of Choptuik 



critical collapse [IC, |11| 



CSS 



DSS 



Considerations and plan. The direct way Figure 2: An illustration of a continuously 
to settle our question would be through nu- self-similar geometry (CSS) as a cone, and a 

1 . 1 , . r ir, J. ij discretely self-similar geometry (DSS) as a wig- 

merical simulation of the system: one would t. j v / t> 

, , . , . , . , , gly cone. The singularity is at the tip. 

need to obtain solutions which are ever closer 

to the merger and with ever higher resolution 

near the high curvature region (where the singularity is about to form). This computation 
would be comparable in difficulty to Choptuik's original discovery of critical collapse |10| 
in that it requires successive mesh refinements over several orders of magnitude. 

However, since we are asking a local question, we may expect or hope that a local 
analysis would suffice, namely the analysis of the metric close to the singularity. On the 
one hand a local analysis is disadvantaged relative to a solution of the whole system in being 
indirect and therefore its results need some interpretation which reduces certainty. On the 
other hand, a local analysis supplies more insight into the mechanism that determines the 
local metric, and it is easier to perform. These latter advantages induced us to prefer the 
local analysis for the current study. The demanding full (numerical) analysis is yet to be 
performed (see however the suggestive but inconclusive results in [p!2| ). 
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Our first assumption is that the local metric is self-similar, as discussed above. Our 
second working assumption is that if several self-similar metrics exist then the one actu- 
ally realized by the system is the metric which is most attractive, or most stable, in an 
appropriately defined manner. 

One local self-similar solution, the "double-cone" has been known for a while Q. In 
terms of the (p, x) coordinates defined in figure || it is given by 

= '^P' + + "^"''^^^ "^^^ + - 4) dnl_^] , (1.1) 



where t denotes Euclidean time (since the solutions are static we may work either with a 
Lorentzian or with a Euclidean signature), and dO|)_3 is the standard metric on the S^~^ 
sphere. Let us recall some of its properties. The (x, t) portion of the metric is (conformal 
to) the two-sphere S^. Thus the metric is a cone over a product of spheres x S^^^ 
which is the origin of its name. Its isometry group is 50(3)^.4 x SO{D — 2)q, which is an 
enhancement relative to the generic S0{2)t x SO{D — 2)q isometry of the system. The 
double cone is smooth everywhere except for the tip p = 0. It is manifestly CSS under 
the transformations p ^ p for any a. Finally, a linear analysis of perturbations around 
the double cone preserving the full SO (3) x SO{D — 2) isometry reveals [|[ an oscillating 
nature (as a function of p) for low enough dimensions, D < 10. 

The current research started from a confusion regarding the CSS/DSS nature of the 
critical merger solution. Q assumed CSS for simplicity and found the double cone. |12| 



found good but not overwhelming evidence for the double cone in numerical solutions. In 
Q a close relation between the merger and critical collapse was discovered,^ raising the 
possibility that just as the critical collapse solution is DSS, the critical merger could be 
DSS as well. Moreover, the linearized oscillations, mentioned in the preceding paragraph 
were realized to be analogous to GHP oscillations (Gundlach-Hod-Piran |l^ , 15 1) in critical 



collapse. In critical collapse the log-period of GHP oscillations is known to be essentially 
the log-period of the DSS critical solution. Therefore Q viewed the oscillations as pointing 
towards a DSS nature. 

Accordingly, our first objective was to find a DSS solution to the system. Since we 
saw little hope in finding an analytic solution, we turned to a numerical method. Rather 
than simulate the whole system and tune a parameter for criticality, we followed [^] and 
imposed periodic boundary conditions. We used two different algorithms. The first was 
of a relaxation type: the 3 fields are solved for iteratively using a selected 3 out of the 5 
equations. This method involves an interesting interplay between the usual local variables 
(the fields) and certain global variables. The essential idea in the second algorithm is to take 
as a merit function the sum of squares of all 5 equations. Altogether, despite considerable 
work the code never converged to a (new) DSS solution, but rather to the double-cone. 
Therefore we relegate the description of algorithm and implementation to appendix ^ and 
choose to detail only the second approach. 



^See also |13| which followed and studied the dimensional dependence of the Choptuik scaling constants. 
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The apparent paradox between the existence of GHP oscillations ad the absence of a 
DSS solution is explained, in hindsight, by the fact that while DSS indeed implies GHP 
oscillations the converse is incorrect: oscillations may arise also around CSS solutions. 

The misguided fruitless search for a DSS solution pointed us towards a different effort, 
which is the focus of this paper: one could analyze the linear stability of the double cone. 
If it is found to be unstable (in a sense to be described below) then it is very unlikely to 
be the metric chosen by the black-hole black-string system. 

Indeed, the asymptotic boundary conditions, including the compact nature of the extra 
dimension, can be viewed as an asymptotic perturbation of the critical merger solution. 
By assumption, this perturbation is irrelevant near the singularity. In addition one could 
consider turning on various perturbations far away from the system, such as putting the 
black object in a non-flat but low-curvature background or turning on a cosmological 
constant (note that these perturbations belong to a wider class - the first does not obey 
the generic isometrics and the second perturbs the equations). If the double-cone is found 
to be unstable to some asymptotic perturbation it would be unlikely to be realized as the 
critical merger solution. On the other hand, if it is found to be stable that would make it 
a viable candidate for being the critical merger solution. 

Stability. The preceding discussion motivates us to formulate our stability criterion. Ac- 
tually, we are seeking a solution which is not absolutely stable to asymptotic perturbations, 
but one which has a single unstable such mode - this is the mode which corresponds to 
motion on the branch of black-string (or black-hole) solutions away from merger.^ There- 
fore we define a self-similar solution to be stable if all but one asymptotic perturbation 
are irrelevant at the singularity, and we proceed to define "irrelevant" and "asymptotic 
perturbation" . 

Each mode can be characterized by its p (see figure |l|) dependence, which must be 
a power law, since the background is continuously self-similar and therefore the scaling 
generator can be diagonalized simultaneously with the small perturbations operator (the 
Lichnerowicz operator) - this will be seen explicitly in section For each mode we define 
a constant s through 



where Sg^i, is the perturbation to the metric, and in general s could be complex. Actually 
since the modes are determined by a system of second order ordinary differential equations 
(ODEs), the modes come in pairs, and we denote the corresponding pair of s constants by 
s±, ordered such that 5R(s_) < 3^(5+). 

We refer to a mode as irrelevant if it is negligible close to the tip (the singularity), 
namely if 



^More precisely, motion onto the two branches, that of a black-hole and that of a black-string, is asso- 
ciated with two modes defined up to multiplication by K+ and these modes are not necessarily related by 
a multiplication by -1. See the last paragraph of section 3 for further details. 




s > . 
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We define the asymptotic perturbations as those corresponding to s+. The rational 
behind the definition is common-place: for example in electro-statics, solutions of the 3d 
Laplace equations come with two possibilities for the radial dependence for each angular 
number I, being either r' or r~'~^. The mode is interpreted as an asymptotic perturba- 
tion, while the r~'~^ is interpreted as a perturbation to the source which lies at the origin. 
Our definition is ambiguous when 5R(si) = 5?(s2), but it will happen only for the / = case 
which is studied in detail in section ^, and will turn out to pose no problem. ^ 

Altogether our definition of stability as the case when all asymptotic perturbations 
but one are irrelevant at the tip means that 



s+ > 



(1.4) 



for all but one perturbation. It can be said that such a solution is a co-dimension 1 
attractor for asymptotic perturbations. We note that this definition differs from other 
common definitions of stability which involve the positivity of the Lichnerowicz operator 
or absence of imaginary frequencies in modes with time dependence. 



Method and results. In section ^ we determine the spectrum of s constants ( |1.2| ) for all 
perturbations with the generic SO {2) x SO{D — 2) isometry. We use an action approach: 
we write down the most general ansatz consistent with the isometry which includes 5 
fields which depend on two variables, compute the action and expand it to second order 
in perturbations. Then we fix a gauge, derive the corresponding pair of constraints from 
the action and the (other) 3 equations of motion. The latter are solved by separation 
of variables and the solutions are then tested against the constraints. Our result for the 
spectrum is summarized in ( 2.31 ). It consists of two families of solutions Sj^{1), s\.{l), a 
scalar and a tensor with respect to S^^^ . 

For all but one^ / = mode we find that S4. > and hence as explained above the 
double cone is a viable candidate to be the critical merger solution. Combining this with 
the fact that even after our search for DSS solutions described in appendix ^ the double- 
cone remains the only known self-similar solution with these symmetries, we interpret the 
result as strong evidence that the double cone is indeed the critical merger solution. 



Non-linear spherical perturbations. The spherical (/ = 0) perturbations are special. 
There are two such modes and in |Q evidence was given that there are two specific linear 
combinations that generates "smoothed cones" (see figure ^) by attempting a Taylor ex- 
pansion of the fields and equations around the smooth tip of an assumed smoothed cone 
and finding no obstruction to a solution. In section |^ we confirm this by an analysis of the 

^While our definition is intuitively clear it difi'ers from the usual definitions requiring smoothness and 
normalizability at the tip. Since the double cone is singular the usual prescriptions do not apply. Concep- 
tually one could determine the perturbation spectrum around the smoothed cone demanding the standard 
smoothness and normalizability at the smooth tip, and then rescale towards the double-cone (see figure^. 
It would be interesting to test whether this limit would result in our "s+ prescription". Another possibility 
would be to perform a non-linear admissability analysis of the modes, which we indeed perform in the I — 
sector, as explained below. 
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qualitative features of the full non-linear dynamical system. The perturbation associated 
with the smoothed cones is precisely the special relevant mode mentioned above that moves 
the solution off criticality and along the solution branch^ . Any other linear combination is 
seen to be highly singular at the tip, which justifies us in discarding it (and it is consistent 
with our "s+ prescription"). 

From the analytical point of view the qualitative analysis shows an interesting feature. 
The system is non-integrable (see the discussion below eq. p.lO| ) . The qualitative dynamics 
of 2d phase spaces is quite limited and never chaotic while in higher dimensions chaos is 
common. The phase space of this system is 3d: there are two phase space dimensions for 
each of the two modes minus a constraint. However, the dynamical system inherits the 
scaling symmetry of the background. One can define a reduced 2d phase space system by 
choosing a plane transversal to the symmetry flow, and then define a reduced dynamical 
flow to be the original flow projected onto the plane through the symmetry flow. This 2d 
phase space system is now amenable to analysis through the determination of equilibrium 
points: focal, nodal or saddle, and we are able to solve for the full qualitative features. The 
solution is summarized in figure ^ 

It turns out that a related dynamical system, given 
by the Hamiltonian H = {p^ + p^) /2 — /2, ^ was al- 
ready analyzed in quite a different physical setting as a 



model for critical phenomena |17] (see also a higher di- 
mensional generalization |l^]). There one analyzes mini- 
mal surfaces in a black hole background and one finds the 
(local) critical solutions to be cones. The appearance of 
essentially the same dynamical system (system of ODEs) 
in different physical settings suggests that this dynami- 
cal system is in some sense the minimal example for the 
physics of criticality including self-similarity and critical 
exponents. Moreover, the same Hamiltonian, apart from 
a sign change of the potential which is inessential for cur- 
rent purposes, was already studied in in the context 
of Yang- Mills theories. 




Figure 3: A single smoothed 
cone solution can be scaled down 
to provide a continuous family of 
metrics which approach the cone. 



2. Perturbations of the double cone 

In this section we compute the spectrum of zero modes around the double cone (preserving 
the isometries of the black-hole black-string system). 



Action. Let us start by considering the following ansatz, which is the most general one 
given the U{\)t x S0{D — 2)q x S2,t isometries of the black-hole black-string system, where 
^2 r stands for time reflection t ^ —t 



fdpf + e^'^dt^ + e 



2C 



*In order to cast this system in a form similar to ours (3.10) we transform first to the Lagrangian 
L = 



3 • 



(2.1) 



(i^ + + a;^ y^) /2 and then change variables into u — Iog(a;), v — \og{y) to obtain L = 



(e^"w^ + e^" i]'^ / v"^ + exp(2u + 2t;)) /2. In the H ^ Q sector we may multiply L (or H) by any lapse 
function, and we choose exp(— u — v) to arrive at L = (e""" v? + e""" w^/u^ + exp(u + v)) /2. This form 



still differs from (3.1C ) but it is similar as it has a nearly canonical kinetic term and an exponential potential. 
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All the fields are functions of p and x only. The ansatz is "most general" in the sense 
that all of the Einstein equations can be recovered by varying the gravitational action 
with respect to these fields. Prom now on, in order to shorten the notation we denote the 
derivative with respect to x by a dot, whereas the derivative with respect to p is denoted 
by a prime. 

After some tedious computation the Lagrangian of the system can be obtained 

S = J exp{^ + Bp + B^) dpdx C 

C = Ki + K2p + K2^-V , (2.2) 

where 



— — i^i := dC d{-^ + ^-C) 



+ (C + - C') + C + $ - C)) / e 



2B, 



e^""" K2n := -2vi/(// + /' + /2R 



-2^' (4 + 2/4) 

+2/(44-44) 

+2/ (24 -4) -2/' (24 -4) 



■■= -2^4 

V := {D- 3)(D - 4) exp(-2 C) , (2.3) 



and we define 



*:=$ + (D-3)C. (2.4) 



The Lagrangian was divided into several parts as follows: V is the potential - a part 
without derivatives and Ki is a kinetic term which is co- variant in the (p, x) plane. The 
rest of the kinetic part was somewhat arbitrarily divided such that terms with a e''^^p 
factor were collected into and the term with a e"^^^ factor was denoted by K2y-. 

Gauge fixing. This action is invariant under reparameterizations of the {p, x) plane (2 
gauge functions). We do not pretend to know an optimal gauge, but we start by fixing 

/ = (2.5) 

as it seems to considerably simplify the equations. The corresponding constraint comes 
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from the equation of motion for / [p, x) and is given by 
1 6S 







2 6f 



-{D -3){<^'C + ^C' - CC) ) e 



-<i> + B,-2B^)e'^~^^+^^ 



-d^[{2B'^-B'y-^^+^^) 



Now substituting the gauge / = into the Lagrangian, we get 



(2.6) 



L 



-2B'^'e 



(L>-3)C"($'-C" + ^')e 



'\^-2Bo 



+{B, 



X ' — ' Bp ' 



{D-?,){D-4)e 



-2C 



^^ + Bn + B^ 



(2.7) 



We choose to fix the remaining gauge freedom (beyond / = 0) such that the kinetic 
term with respect to x i^i^x^ terms) in is canonical (or more precisely, field inde- 
pendent) 

B^ = ^ + Bp + h . (2.8) 
where h = h{p, x) will be fixed later. The associated constraint is given by 
5S 







dBy 



fix Bt, 



2Bp^ + {D-'i)C{^ + ^-C) 
2{^' + B'p + h')^' + {D- 3)C" ($' - C" + ^' 



(2. 



The background. The double-cone is a Ricci-flat metric given by^ 

^2 



P 



D-2 



{dx^ + cos^ X dt^ + {D-'i) dnlo- 



(2.10) 



As mentioned in the introduction, our objective here is to compute the spectrum of per- 
turbations around this metric, namely to solve the Linearized Einstein equations around 
this background. 

The gauge choice ( ^^J2.8|) requires to redefine the angle X ^ X- Choosing 



h{p) = -{D-3)ln l^p 

implies the following p-independent equation for x 

1 1 

m — = m cos X 



D-A 
D-2 



(2.11) 



(2.12) 



^See also (1.1) 
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The solution of this differential equation is given by 

X = arctan(sinh x) (2-13) 

where for simplicity we chose the constant of integration to be zero. Note that while 
X ranges over [— 7r/2,7r/2] x ranges over [—00, +00]. With this redefinition at hand the 
solution ( p.lOP becomes 

is' = V + 1/!^ + - 4) <in|„-,) , (2.14) 

here for simplicity of notation we omit tilde above x ^-^id denote x by X- 



Linearized equations. Now we slightly perturb the double-cone solution while keeping 
the gauge condition (|2.8|) unchanged, that is we set 



5. 



+ 6, = + 6, 



P 



^ ^ \coshx/D^ 



In 



P 



cosh xV-D — 2 



+ 



C = C('^) + c = ln(,,/§^ +c 



(2.15) 



together with the gauge condition derived from (|2.8D 



p ' 



where 



:= {D-2,)c + . 



(2.16) 
(2.17) 



(see also (|2j)). 

Let us substitute equations ( 2.15|J2.1^ ) into the /-constraint ( |2.6| ) and 6-^-constraint 
( |2.9| ). The zeroth order in small perturbations vanishes in both cases, whereas the first 
order terms satisfy the following equations 



{D - 2)bp - pij' - ptanhx (V'' -</>' + = 



(2.18) 



-2p^^" + 2p 



{D-2)b'p-{D-l)i;' +{D-2) 
+2 {D -2){D- 3) {hp - c) = 



^'p + V' 



sinh (2x) 



(2.19) 



Currently we are in the position to derive the equations of motion for the rest of the 



fields. As a first step we substitute ( 2.15 , 27Lq ) into the Lagrangian (2^), expand it and 
keep the quadratic part in the perturbations 
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cosh X 



^-3 



{D-2) 



{3D - 10) {D-3)c ^ + 6(/)V' - 40' ^ + 2bp ' V' 



(L> - 3)c{ip - c + 0) + 26pV' cosh^ X - 4p (2^-' + 6;;,) V 
-2{D - 1)^2 _ 2(D _ 3) _ c + 6p)2 



(2.20) 



As a result, the equations of motion for the fields bp, c and cp are given respectively by 



2 [{D - 3)bp - 2v + = j§:r^p^' + i^ij^" + ^ ^ 

p [(I) - 1) (3^-' - c') - (Z) - 3)6^ + {3r - c" + b'l,) 
+ {D-2) cosh^ X (^V^ - c + 6p) 

2 (D-2) (c-2'iP + (t) + {D -A)bp 



Q = p[{D- l){3^' - 4>') -{D- 3)b'p] + p" {3r - 4>" + b"p) 

+ {D-2) cosh^ xU-'<i> + bp) -'^{D -2){D -3) {bp - c) 



(2.21) 



Let us denote the constraint ( 2.1q ) coming from the gauge fixing condition / = by 
A, and by B the constraint ( 2.19 ) coming from the gauge condition ( ^.161) . One finds that 
the following relations hold on the solutions of equations of motion ( |2.21j ) 

pdpB {p, x) 



^-^^P^^^ = - 2(Z?-2)cosh2x 
d^B {p, x) = 2{D-2)A {p, x) + 2pdpA (p, x) + 2 tanh x B {p, x) 



(2.22) 



It means, that once the constraints are satisfied for some value of the coordinate x they 
will vanish identically for any x- 

2.1 Solving the equations 



We proceed to solve first the equations of motion (|2.21| ) , and then later select those solutions 



which satisfy also the constraints ( 2. IS , 2. 19 ). 



In order to solve the equations of motion ( p.21|) we attempt to separate the variables 
through the following general ansatz 

/b^y 



X{p,x) 



0(1) =XV^P«(tanhx) 



(2.23) 



,(1) 



where I is any non-negative integer, are unknown constant vectors and Pi (tanhx) are 
the Legendre polynomials which can be shown to satisfy the following differential equation 



cosh^ X -r^Pi (tanh x) + / (/ + 1) Pi (tanh x) = 
dx^ 



(2.24) 
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The description "general ansatz" needs some explanation. Once we add a sum over I to 
the r.h.s of ( 2.23| ) this becomes the most general decomposition since and Pi form a 
complete basis of functions. The question is whether the variables separate in ( 2.21J ), thus 
allowing to omit the sum. A priori the 5*0 (3) isometry of the background tells us that 
the angular coordinates can be separated into spherical harmonic functions. In general 
such functions are labelled by l,m but the U{l)t isometry implies that only m = terms 
contribute. If our fields included only scalars then the equations would be guaranteed to 
separate into the scalar spherical harmonics Yiq, namely the Legendre polynomials. In our 
case the perturbation includes also tensor modes (with respect to the S^), but still direct 
inspection ^ confirms that equations ( |2.21| ) separate under the ansatz ([2.231 ). 
As a result we get the following set of algebraic equations 



[s (s -1)E2 + sEi -l{l + l)K]X = VX 
where we have defined the following 3x3 constant matrices 

-{D-3) 2{D-1) 3{D-3){D-1) 
{3D - 5) {3D -5){D- 3) 

-{D-3) 3{D-1) (D - 1) (3D - 10) 

1 2 3(^-3) 
1 {D-3) 
13 {3D - 10) 



(2.25) 



Eo 





1 


{D- 


3) 


K = {D-2) 


1 


{D- 


3) 




1 1 


{D- 


4). 




{D 


-3) 





V = 2{D-2) 


(D 


-3) 


- 1 




AD 


-4) 


- 1 



(2.26) 



-{D-3) 
-2{D-3) 
(7 - 2D) 

The spectrum of s is determined from the following characteristic equation 

Det [s {s-l)E2 + sEi-l{l + l)K -V]={) (2.27) 

After some tedious algebraic rearrangement, one can simplify the above equation and get 

+ (L» - 2) s - (/ + 2) (/ + 1) (D - 2) = 
+ {D-2)s-{l + 2){l-l){D-2) = {) 
s^+CP - 2) s - / (/ - 1) (D - 2) = 



X appears in (2.21) only through the combination cosh x and after use of (2.24) all x dependence 



disappears. 
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The solutions are given by^ 



si,2 = -\2-D± yJ{D - 2) (4/2 + 12/ + L» + 6) 



S3,4 



i f 2 - L» ± ^/{D - 2) (4/2 + 4/ + D - 10) 



i f 2 - D ± V(L» - 2) (4/2 _ 4/ + D - 2) 



whereas the corresponding eigenvectors are 



X 



1,2 



(/ + 2) 4 - 1) ± ^{D - 2) (4/2 + 12/ + D + 6) 



2/ - D + 6 =F ^(1) - 2) (4/2 + 12/ + D + 6) 
2(/ + 2) 



(2.28) 



X 



3,4 



4(L»-2)(L>-3)(/ -!)(/ + 2) 



(Z) - 3) Z) - 2 ± v^(L> - 2) (4/2 + 4/ + D - 10) 
-4(i?-2)(/-l)(/ + 2) 



(/ - 1) A-D± ^{D - 2) (4/2 -Al + D-2) 
2l + D-4± y/{D-2) (4/2 - 4/ + D - 2) 
2(/-l) 



(2.29) 



It turns out that X^^^ satisfy the constraint equations (2.1£,2.18) as well, and thus 
they are part of the perturbation spectrum. On the other hand, Xi^2 and X^^q do not 
satisfy the constraints independently. Since mixing of different modes is allowed for those 
values of s which are degenerate, one concludes that in order to find out other possible 
solutions of the perturbation spectrum we need to take a superposition of Xi^2 and X^^ 
corresponding to the same values of s and then check whether such a combination can 
satisfy the constraints. 

According to ( 2.28| ) one can see, that in order to make the powers of p equal we should 
consider the following superpositions 



1^2 



X'+^Pl+2 (tanhx) + FX{Pi (tanhx) 
X'^^Pl+2 (tanhx) + GX'2Pi (tanh^) 



(2.30) 



where Fand G are constants to be determined. Substituting these linear superpositions into 
constraints equations, we find that the constraint can be satisfied by taking F = G = —1, 
and we have another family of solutions. 

^Indices 1,3 and 5 correspond to upper sign, whereas indices 2,4 and 6 to lower one. 
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In summary, the full perturbation spectrum is given by 



4 = H2 - -D ± y^{D-2) {4P + 121 + D + 6) 
s± = H 2 - ± V(^-2) (4/2 + 4/ + D - 10) 



(2.31) 



where^ / > 0. 

For the modes are given by 



X'+'Pl+2 (tanhx) - X[Pi (tanhx) 
Xl (p,x) = /- \x'+^Pi+2 (tanhx) " (tanhx) 



while for they are given by 



X% {p,x) = X' /j^iP^tanhx) 



(2.32) 



(2.33) 



According to the "s+ prescription" boundary condition (below (|1.3D ), we should elim- 
inate all the S- modes. For s*(/ = 0) this prescription is ambiguous, but after studying 
the I = sector in detail in the next section we will conclude that still the b.c. reduce the 
dimension of the solution space from 2 to 1. 

In addition, another mode should be eliminated from the above perturbation spec- 
trum, namely {p-,x) foi' s+(^ = 1) = 0) since it corresponds to a residual gauge of an 
infinitesimal shift in x coordinate. 

Altogether it can be seen that except for s^(/ = 0) all other s+ (physical) modes are 
positive, thus satisfying our stability criterion ([1.41) . 

3. Non-linear spherical perturbations 

In this section we obtain the qualitative features of the dynamics of the full non-linear 
perturbations in the "spherical" sector (/ = 0, namely preserving all isometries). 



Action. The most general D-dimensional metric with SO{m + 1) x SO{n + 1) isometry, 
D = m + n-|-lis 

ds' = e^''''dp' + e^^dnl + e^^dnl, (3.1) 

where all three functions Bp, A and C depend on p only, dO,"^ and dO,"^ are the standard 
metrics on the m and n spheres, and there is a reparameterization gauge freedom p p'{p)- 
For applications to the black- hole black-string system we need only the case m = 2, which 
represents the Xi t 2-sphere while the n-sphere is the angular sphere. 
The action is 



S 
K 

V 



dpe 
mn 



m A+n C ^Bp 



V 



—m{m 
m (m - 



-{A'-C'f-(l-^^ {mA' + nC) 
1)C7'2 - 2mnA' C 



DA 



n[n 



l)e"^^ + n(n 



l)e 



-2C 



t" stands for "tensor" and "s" for "scalar". 



(3.2) 
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where a prime denotes a derivative with respect to p (the overall sign was chosen such that 
S = -j^R). 

The system enjoys a scaling symmetry ^ 

ds^ e^'^ds^ (3.3) 

namely 

Bp-^ Bp + a 
A + a 

C ^ C + a . (3.4) 

It is convenient to fix the gauge such that the kinetic term is canonical (more precisely, 
its prefactor is field independent), namely 

Bp = mA + nC . (3.5) 

The ansatz reads 

ds'' = e2-^+2nC^^2 ^ e'^^dnl + e^^d^l , (3.6) 
while the action becomes 

S = J dp[K-V] (3.7) 
V := exp(2mA + 2nC)y = m(m-l)e2('"-^)^+2"^ + n(n-l)e2"*^+2("-^)^ , 
and it is supplemented by the constraint 

= H ■.= K + V . (3.8) 

Change of variables. It is convenient to make the following field re-definitions. First one 
makes a linear re-definition that simplifies the potential term 

u := 2(m - 1) ^ -I- 2n C -I- log (m(m - 1)) 

v:= 2mA + 2{n-l)C + log (n(n - 1)) . (3.9) 
The potential and kinetic terms become 

= e" + 

K 



A{D-1) 
1 

A{D - 2) 



mn {u' — v')'^ — — — -r {mu' + n v')'^ 
V " ~ 2) 

[m(n — 1) u"^ + n{m — 1) v''^ — 2mn u' v'^^ (3.10) 



^The variation of the action due to an infinitesimal symmetry operation is S^S — S. Usually one 
considers symmetries which vary the action only by a boundary term thus defining a conserved current, 
which is absent in this case. Still this variation is enough to guarantee that the equations of motion are 
satisfied. Indeed, assuming we have a solution for the equations of motion 5S/S<f> — (where is a collective 
notation for the fields) then their variation vanishes as well 6a SS/d(f> = 5{daS)/S4> = SS/Stj) = 0. 
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This is a system with two degrees of freedom and a potential which is a sum of two 
exponentials. At first it appears to be similar to a Toda system, but actually it is probably 
not integrable for the following reason: the spring potential in a Toda system is of the form 
— X (where x measures the deviation from equilibrium length) , and while the linear term 
often cancels due to the mass being acted on by two springs, one from each side, here there 
are only two springs for two masses, and thus our system which lacks a linear term is not 
of a Toda form. 

Still it is convenient to make the following "Toda inspired" change of variables 



X2 



Pi :-- 
P2 :-- 



u 



(3.11) 



where the notation should not be mistaken to imply that the X^s and P's are conjugate 
variables. The equations of motion read 



X[ 
X'2 



PiXi 
P2X2 

and the constraint ( |3.8D becomes 
1 



P{ = 2((l-i)Xi+X2) 
P^=2(Xi + (l-i)X2) ' 





The {Xi,Pi] 



4{D-1) 



mn{Pi- P2)' 



{mPi+nP2f 



+ X1+X2. 



(3.12) 



(3.13) 



1, 2 variables also provide a convenient realization of the scaling symmetry 



Xi - 




Pi - 


^ e'^Pi 






p - 





(3.14) 

where a is related to a in ( |3.4[) through a := {D — 2) a. 

The double cone metric in the X, P variables is found either by transforming ( |1.1| ) 
according to the changes of variables (^, 3.11 ) or by solving directly the equations of 
motion ( pl|) subject to scahng (|1|) invar iance. It is given by 

2 



Pi = P2 



m 



D 



n 



D 



P 
1 

27 
1 

27 



(3.15) 



Next we transform to 



X+ :-- 

X- :-- 

P+ : = 

P :-- 



X1+X2 
1 



Xi 



1 



X2 
m n 

mPi + nP2 
P1-P2. 



(3.16) 
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To arrive at these definitions we first considered tfie 5d case m = n = 2 where by symmetry 
it is useful to transform to Xi it X2, Pi =t P2 and then we generalized to arbitrary m,n 
paying attention to the form of the kinetic energy ( 3.10| ). The equations of motion read 

X+ ' = — 3— (X+ P+ + mn X- P-) 
m + n 

X~ ' = — — (X+ P^ + X- P+ + (n - m) X- P') 
m + n ' 

P+' = 2{D - 2)X+ 

p-' = -2X- (3.17) 



Fixing the scaling symmetry. Now comes the crucial step in the analysis, which will 
allow the qualitative solution of the dynamical system. The phase space consists of 4 
variables Xj, Pi constrained by (|3.13| ) and hence it is 3d. However, dynamical systems in 
3d can be quite involved, and we would not know how to analyze this system. Fortunately 
the scaling symmetry (^, 3.14 ) can be used to reduce the problem to a 2d phase space, 



where the number of qualitative possibilities is quite limited and a full qualitative analysis 
is possible. 

The idea is to fix the symmetry by choosing a 2d cross section of the phase space which 
is transverse to the symmetry orbits. Then we supplement the infinitesimal p-evolution by 
an infinitesimal symmetry operation such that we always remain on the 2d cross-section, 
thereby reducing the problem to a 2d phase space. 

In practice we fix the symmetry as follows. Being transverse to a scaling symmetry 
means introducing an arbitrary scale. We choose the following condition 

X+ = 1 . (3.18) 

In order to define the reduced "scaling compensated" evolution we introduce a con- 
densed notation for the phase space variables 

r = (Xi,X2,Pi,P2) . (3.19) 

The scaling transformation is given by y* e'^'" y* where 

g' = (2,2,l,l) . (3.20) 

We denote the functions on the r.h.s of ( 3.12| ) by P{Y^) such that 



Y" = f{Y^). (3.21) 

Note that since p carries scaling charge qp = —1 /* carry charge + 1. 

We need to distinguish the reduced evolution parameter pji from p since scaling acts 
on p as well dll). Now we can define the reduced phase space trajectory = Y^{pr) 
by 

Yk = r-<f4^Y^ (3.22) 



dpR " ^ 2X+ 
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where f~^/{2X~^) is the compensating infinitesimal scahng parameter and is defined such 
X'^ = 0. For clarity, we write down the definition of explicitly 

f^^^R^ -■= +^^^rPr)- (3-23) 

We parameterize the reduced 2d phase space by {X, 9) where 

X^X],, (3.24) 

and is a hyperbolic angle which parameterizes the hyperbola in P space given by the 
constraint ( 3.13| ) together with the condition ( |3.18| ), namely 



cosh e = ^ = P+ 

2^{D-2) {D-l) ^ 



^inh^ = -^V^ Pr ■ (3.25) 



The signs above depend on conventions: the sign in the definition of cosh(^) is related to 
setting the direction of the flow towards the tip, and the sign in the definition of sinh(0) 
sets the sign of 6. Our (final) expression for the reduced equations is 

d 2 

X = j= sinh6l(l + nX) (1 - mX) 



dpR ^mn 
d 



dpR 



6* = VI)^sinh6' + cosh6' X , (3.26) 



where we chose to rescale pR Pnl^rn + n. 

Given a solution of ( 3.26 ), or equivalently functions Y^{pji) satisfying ( 3.22| ), we still 



need to uplift it to a solution of ( 3.21| ). First we integrate for the scale factor evolution 
a{pR) 

Then we define the uplifted trajectory Y^{p) by 

dp = e~" dpR . (3.28) 
A direct computation confirms that with this definition ( |3.21| ) is satisfied 

dp dpR dpR dpR 



where in the first equality we used (|3.28| ), in the second to last we used (^ , |337| ) and 



finally in the last equality we used the fact that /* has charge + 1 with respect to scaling. 
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3.1 Analysis of the reduced 2d phase space 

The analysis of the quahtative form of the reduced 2d phase space ( f3.26|) proceeds by 
determining the domain, equihbrium points, their nature (at hnear order) and a separate 
analysis of the behavior at infinity. 

Domain. From the definition of Xi, X2 ( p.ll|) it is evident that they are both positive. 
Given the gauge fixing condition (3.18) and the change of variables ( 3.16| ) we find that the 



domain is restricted to the strip 

-- <X < — . (3.29) 

n m 

The boundary is strictly speaking not part of the domain, but since the equations ( |3.2(^ ) 
continue smoothly to the boundary while X vanishes there, we can join it to the domain 
and allow for equalities in ( p. 29 ). 



Equilibrium points. There are 3 finite equilibrium points 

Xn = 6(^ = 



X, = ± tanh(0i) = -V^^ , (3.30) 
X, = -i tanh(02) = ,M ^ 



two of which are on the boundary. The point (0, 0) corresponds to the double cone ( 3.15| ), 



which is indeed a fixed point of the scaling symmetry. The role of the other points will 
become clear below. 



Linearized analysis of equilibrium points. The linearized dynamics around the an 
equilibrium point {Xq,9q) is given by the 2x2 matrix L 



d 


6X 


= L 


5X 








dpR 


59 




56 



(3.31) 



where 5X := X —Xq, 59 := 9 — 9q. We remind the reader of the classification of equilibrium 
points 

• If the eigenvalues of L are complex (and necessarily self-conjugate since L is real) 
then trajectories circle {Xq,9q) or spiral around it, and the point is called a focal 
point. 

• If the eigenvalues of L are real and of opposite sign then {Xq,9q) is called a saddle 
point and it is unstable in the sense that only if the initial conditions are fine-tuned 
to lie on a specific trajectory will the evolution flow into this point. 

• If the eigenvalues of L are real and of the same sign then {Xq,9q) is called a nodal 
point. If they are both positive then the point is called repulsive, while for negative 
it is called attractive. Naturally, the adjectives repulsive and attractive interchange 
under inversion of the flow (time reversal). A focal point may be called repulsive (or 
attractive) as well if the real parts of L's eigenvalues are all positive (or negative). 
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-0.40.2 0.20.4 



Figure 4: The phase space (for m = n = 2). There is one interior equihbrium point at (0,0) 
which represents the double cone. For _D < 10 it is focal repulsive but the spirals cannot be seen in 
the figure since their log-period is too large, while for Z? > 10 it is nodal repulsive, and hence the 
figure represents the whole range of dimensions. There are two finite equilibria on the boundary at 
(I/to, 0i) and (— 1/71,^2) which are saddle points. Each one has a critical curve which approaches 
it (heavy line). These two special trajectories denote the smoothed cones, where either S" or 
gm g]^pjj^]jg smoothly while the other one stays finite. Finally there are two attractive equilibria at 
infinity at (l/m, —00) and (— 1/rt, -I-ck)), and the thin lines trajectories represent generic trajectories 
which end at these infinite equilibrium points. 



For the double-cone {X, 9) = (0,0) we find 







(3.32) 
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The eigenvalues of this matrix are 



A± = - ( VD-2 ± VD - 10 



(3.33) 



Let us make several observations 



The eigenvalues depend only on D and not on m, n separately. 
= 10 is a critical dimension, as found in 

— For D < 10 the eigenvalues are complex and we have a repulsive focal point. 

— For D > 10 the eigenvalues are real and we have a repulsive nodal point. 

— For the critical value D = 10 the eigenvalues degenerate, but L is not pro- 
portional to the identity matrix, but rather L can be brought by a similarity 
transformation to the form 

^ y/2 1 
V2 

• In the D ^ oc limit we have A+ ~ \/^, A_ ~ 2/^/D. 

• The eigenvectors which correspond to the eigenvalues X± are [2, — X± \Jmn\ 
At the equilibrium point {X\ ,9i) we find 





(3.34) 



Li = cosh 01 



-2{D-1) 

\/mn 



(3.35) 



Similarly L2 is obtained by substituting 9i ^02, ?7i <-> n. Since detL < we see that at 
both {l/m,0i) and (— 1/n, ^2) we find a saddle point. The positive direction (eigenvector 
of the positive eigenvalue, defining the repulsive direction) is [0, 1] in both cases, namely 
along the boundary. 

Behavior at infinity. There are two attractive equilibrium points at infinity 

X'^ = — 0'^ 



— 00 
+ 00 



(3.36) 



(^3,63) is attractive since for < 0i we have dX/dpji > 0, d0/dpii < 0. There are no 
other fixed points at = —00 but with —1/n < X < 1/m since liuig^^oo dX/d0 < 0. 
Moreover, close to (^3,^3) the trajectories are found to asymptote exponentially fast to 
X3, namely \5X\ ~ exp(— const 0). Analogous properties hold for the attractive equilibrium 
point (X4, 04). 

The whole picture. Having found all the ingredients above we may assemble them into 
a complete phase diagram, figure which is discussed in its caption. 

The trajectories which end at {Xi,0i) and (^2,^2) represent the smoothed cones. 
This is seen by transforming the definition of the smoothed cone through the changes of 
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variables, as we proceed to explain. Before gauge fixing, one of the smoothed cones can be 
defined to behave in the limit p — > as i?p = 0, A ~ log(p), C ~ Co (while for the other we 
interchange A <^ C). After changing gauge and transforming to the {X,P) variables this 
becomes Xi = m/(m— X2 = const p~'^"^/^"^~^\ Pi = —2/p, P2 = —2m/{m—l)p^'^ 
which in the {X,6) variables tends to (Xi,9i) as defined in ( |3.30| ). 

In summary, we have confirmed non-perturbatively the existence of the smoothed 
cones, going beyond the perturbative analysis of |^] near the smoothed tip. Altogether 
there is a one-parameter family of solutions which in the linear regime (far away from 
the tip) corresponds to the various linear combinations of the two linearized modes. Two 
of these solutions correspond to the smoothed cones, while all the rest end at infinity in 
phase space and yield a singular geometry. Even though strictly speaking our b.c. "5+ 
prescription" is undefined for / = since 3^(si) = 5R(s2), we define it to mean to retain 
only the smoothed cone solutions, thereby the double-cone is an attr actor at co-dim 1, as 
wanted, rather than at co-dim 2. Equivalently, these two special directions are defined 
up to multiplication by M+, and hence we loosely refer to them as a single mode which is 
normally defined up to multiplication by M. 
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A. A numerical search for a DSS solution 

In this section we present the details of our frustrated search for a Discretely Self-Similar 
(DSS) solution to the system, which if found, would have been a candidate to be the critical 
merger solution. 

Formulation of the problem. The Lagrangian for the metric in the following form 

ds'' = e^^'^^P'^^dp^ + e^^'-^P'^^dx - f{p,x)dpf + e^''^"^^^ dt" + e^<P^^U^l_^ , (A.l) 
where p here is the log of p in ( |2.1[) , and choosing the gauge 

f{p,x) = 

Bp{p,x) = B^{p,x)^B{p,x) (A.2) 

is 

L = -e* [2B'^' + {D- 'i)C'{^' + - C) + (' ^ )] - 

(L»-3)(L»-4)e*+2^-2^, (A.3) 
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where ^' = $ + (Z) — 3)C as in ( |2.4D and prime and dot denote derivatives w.r.t p and x- 
The constraints obtained from the gauge fixing ( [A. 2D are 



^' + {D- 2,)C' + + {D- 2,)C'C - (B'^ + B<^>') - {D - 3){B'C + BC') = 0, (A.4) 

- $ + (D - 3){C" -C) + - $2 + (i:) _ 3)((7'2 _ ^2)_ 

- 2(5'$' - B^) - 2{D - 3){B'C' - BC) = 0, (A.5) 

and the equations of motion which follow from the Lagrangian are 
V^^ = -{D - 3){C'^' + C^) - ($'2 + 

V^C = {D- 4)e2^-2c _ - 3)(C"2 + C^) - (C'^.' + C$) (A.6) 

where is a 2-dimensional Laplacian = dpp + d^^^. 

We are looking for a Discretely Self-Similar (DSS) solution to these equations, one 
which satisfies g^uip + ^piX) = 9iJ.v{p-,x) some A and Ap. In our equations we 
change B ^ B + -^p, ^ ^ ^ + -^p, C ^ C + -^p. The fields thus defined are periodic 
with period Ap. Finally we make one more redefinition <^ — > $ + lncosx/2, influenced by 
the form of the double-cone solution (|1.1|) . With the flelds thus defined the double cone 
solution is constant 

C(p.x) = ^lnff|. (A.7) 

The boundary conditions which the fields should satisfy are 

1. Reflection symmetry at the x = axis 

5(p,0) = $(p,0) =C(p,0) =0 (A.8) 

2. Regularity at the horizon which is at x = There are 5 regularity conditions (3 
from EOMs and 2 from constraints), but only 4 of them are independent. In terms 
of the newly deflned flelds these conditions read 

= 5(p,7r) = ci>(p,7r) =C(/>,7r) , 

S'(p,7r) =$'(/,, vr) . (A.9) 

3. Periodicity in the p direction as mentioned above 

= $(p + Ap) - $(p) = B{p + Ap) - B{p) = C{p + Ap) - C{p) (A.IO) 
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Since our equations are invariant under the symmetry of adding 2 independent constants 
to the fields: B ^ B + consti, C ^ C + consti, $ ^ <I> + const2, the last regularity 
condition can be integrated to the form B{p,7r) = <I>(p, vr) where setting the integration 
constant preserves only one remaining symmetry (scaling) 

B ^ B + const, C ^ C + const, ^> ^ ^> + const. (A.ll) 



Algorithm. One way to find a DSS solution numerically is the following. Using the 5 
equations that we have (3 elliptic and 2 hyperbolic) written as 



fi{B{p,x)MP,x),C{p,x),^/^p)=^ 
with i = 1..5, one constructs a "functional" 



(A.12) 



Fo\B{p,x)MP,x),C{p,x),^,\]--=iZ I'^^l dpf!{B{p,x),Hp,x),C{p,x)). (A.13) 



Note that in addition to the 3 (local) fields, the functional depends also on 2 global variables 
A and Ap. This functional vanishes if it is evaluated on a solution and otherwise it is 
positive, so one can minimize it numerically. We are interested in a solution that satisfies 
our boundary conditions. Those in the x direction (A.8, A.9| ) can be incorporated in the 
algorithm by defining one more non-negative functional: 



Fi[B{p,x),HP,x),C{p,x)] ■■= J dp B\p,0) + ^\p,0)+C\p,0) 



Ap 



+ 







B\p, vr) + <l>\p, vr) + C\p, vr) + {B{p, vr) - $(p, vr))^ . (A.14) 



The boundary conditions in p direction ( A. 10 ) is incorporated differently, as will be ex- 
plained below. 

The remaining symmetry ( [A.ll| ) in the equations can be fixed in many possible ways; 
the most symmetric one being to set the overall average of the three fields to 0. It is done 
by the following functional: 



F2[B{p, x), Hp, x), Cip, x)] ■=[ dx dp(B{p, x) + ^(p, x) + C{p, x))] ■ (A.15) 





The total functional to be minimized is 



F :— Fq + Fi + F2 



(A.16) 
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In order to prevent the procedure from converging to the double cone which is charac- 
terized by constant fields, we multiply F by \{vp + ^) where the "variability" Vp is defined 

by 

TT Ap 

Vp ■■= j dxj dp[B'\p, x) + ^'\p, X) + C'\p, x)) . (A.17) 



The Vp term is added to prevent an endless drift towards an opposite limit, that of rapidly 
varying fields (since if f p ^ oo then ^ would go to no matter what happens to F) and a 
factor 2 is added in order to fix the normalization since the value of the function x + 1/x 
at its minimum is 2. 

The arbitrariness p ^ p + const is harmless, as will be explained below. 

Discretization. By discretizing the domain F becomes a function of a large number of 
(local) variables: F = F(Bij,^ij,Cij) where Bij, ^ij, Cij are the values of fields at the grid 
sites. We make the grid cells "non-square" as follows: if the step-size in x direction is h 
then we take the step-size in p direction to be ^ • c? with d being adjustable parameter (in 
fact if the number of grid sites in both directions are equal then d = ^fi)- In order to 
impose the boundary condition in the p direction we calculate derivatives w.r.t. p at the 
boundaries of the grid with values of fields at the opposite side of the grid, thus making 
our domain effectively cylindrical. 

One can start with some configuration of fields and change them together with the 
global variables A and Ap so that F is minimized (remembering that only F = minima 
arc solutions). This is done by calculating the gradient dF/dBij, dF/d^ij, dF/dCij, 
dF/dA, dF/dAp (the direction in a "space of parameters" in which F grows in the fastest 
way) and going a bit in the opposite direction. This procedure is repeated until we arrive 
at a minimum. There is an additional parameter which sets the "speed of relaxation" 
by defining the change in any parameter atobea— *-a — J7^. Instead of working with 
the global variables A and Ap we find it convenient to use the equivalent variables d 
and 5 = Finally, the symmetry p ^ p + const is of no importance because nothing 
in the algorithm produces running in this direction. Moreover, discretization makes this 
symmetry discrete and creates a finite distance between nearby degenerate vacua. 

Results. We ran a program with D = 5, and so far the only solution found is the 
double cone, even after using various "tricks" designed to avoid it, as described below. If 
multiplication by ^(?;p + ^) is not used the program converges to the double cone very fast 
no matter what the initial configuration is. If this factor is included then it still converges 
to this solution with d — > oo. If one adds to F an additional factor \{d + ^) then it 
still converges to the double cone but in a much longer time, so that F tends to faster 
than l/vp tends to infinity (the configuration becomes very sensitive to changes, so it is 
necessary to take very small, for a grid lOx 10 it should be of order 10~^ whereas without 
this additional factor it works well with $7 k. 10~^). In all cases 5 becomes close to that of 
the double cone solution k, 0.2887. 
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